Magnetic field screening and mirroring in graphene 
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The orbital magnetism in spatially varying magnetic fields is studied in monolayer graphene within 
the effective mass approximation. We find that, unlike the conventional two-dimensional electron 
system, graphene with small Fermi wave number &f works as a magnetic shield where the field 
produced by a magnetic object placed above graphene is always screened by a constant factor on 
the other side of graphene. The object is repelled by a diamagnetic force from the graphene, as if 
there exists its mirror image with a reduced amplitude on the other side of graphene. The magnitude 
of the force is much greater than that of conventional two-dimensional system. The effect disappears 
with the increase of &f- 



Graphene, an atomic sheet of graphite, has a pecu- 
liar electronic structure analog to a relativistic parti- 
cle, and its unique properties have been of great in- 
terest. In graphene the conduction and valence bands 
stick together with a linear dispersion. The low-energy 
physics is successfully described by the effective-mass 
Hamiltonian analogous to the massless Dirac Fermion, 
0, [3, & |3 and its unique transport properties were stud- 
ied. [Jjjp, d, 0, H] Since its experimental fabrication, 
@, IToL llll| the graphene and related materials have at- 
tracted much attention and have been extensively inves- 
tigated in experiments and theories. 

The electronic property of graphene in a magnetic field 
was first investi gate d in theories as a simple model of the 
bulk graphite. [12j There it was shown that, in a uni- 
form magnetic field, graphene exhibits a huge diamag- 
netic susceptibility due to the orbital motion of electrons, 
which is quite different from the conventional Landau 
diamagnetism. The orbital magnetism was also studied 



layered graphenes, [23 



for related materials such as the bulk g raphite . [13l . Il4 1 
graphite intercalation compounds, [la, I la . 1 1 711 ca rbon 
nanotube, 18, 1|| disordered graphene, [2 ' 
" ~ 



I2H, 122J few- 
and nodal Fermions [26j |. 
Quite recently, the graphene in a spatially modulated 
magnetic field was studied in the context of the electron 
confinement using magnetic barrier. [27], [28|, [2i| The dia- 
magnetic susceptibility was experimentally observed for 
quasi-two-dimensional graphite, which is a random stack 
of graphene sheets. [3Q] 

In this paper, we study the orbital diamagnetism in 
non-uniform magnetic fields in monolayer graphene. Us- 
ing the effective mass approximation and the perturba- 
tion theory, we calculate the electric current induced by 
an external magnetic field with wavenumber q, to obtain 
susceptibility x(q) for general Fermi energies. We apply 
the result to arbitrary geometries where a certain mag- 
netic object is located above graphene, and estimate the 
response current induced on graphene, as well as the dia- 
magnetic repulsive force which works between graphene 
and the magnetic object. We find that graphene has a pe- 
culiar property of magnetic mirroring, where the counter 
field induced by the response current mimics a mirror 



image of the original object. 

We start with the general formulation of the electric 
response to the spatially-varying magnetic field in a two- 
dimensional (2D) system. We assume a uniform 2D sys- 
tem on the xy plane, and apply a magnetic field per- 
pendicular to the layer B(r) = [V x A(r)] z with vector 
potential A(r). Here r = (x,y) denotes 2D position on 
the graphene while we later use p = (x,y,z) to spec- 
ify the point in three-dimensional (3D) space. We define 
j(r) = (j x ,jy) as the 2D electric current density induced 
by the magnetic field. Within the linear response, the 
Fourier-transforms of j(r) and A(r) are related by 



(i) 



with response function K^. The gauge invariance for A 
requires J2 V i-T M „(q) q v = 0. The continuous equation in 
the static system, V-j(r) = 0, imposes another constraint 
^2 Kftv^q) = 0. To meet both requirements, tensor 
needs to be in the form, 



(2) 



On the other hand, because V-j(r) = 0, we can express 
j(r) as j x — cdm/dy, j y = —cdm/dx, with m(r) being 
the local magnetic moment perpendicular to the layer, 
and the light velocity c. In the linear response, its Fourier 
transform is written as 



m(q) = x(q)-B(q), 



(3) 



with the magnetic susceptibility x(q). Equations (Q]) and 
([3]) are complementary, and both response functions x(q) 
and K(q) are related by 



X(q) = -A^(q) 

cq z 



(4) 



Graphene is composed of a honeycomb network of car- 
bon atoms, where a unit cell contains a pair of sublat- 
tices, denoted by A and B. The conduction and valence 
bands touch at the Brillouin zone corners called K and 
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K' points, where the Fermi energy lies. The effective- 
mass Hamiltonian near a K point in the absence of a 
magnetic field is given by 0, 0, 0, 0] 



Ho = Hv 











hvk ■ a, 



(5) 



where v is the constant velocity, k = (k x , k y ) = — iV, and 
<t = ((Ta-jCTj,) are the Pauli matrices. The Hamiltonian 
([5]) operates on a two-component wave function (ipA, iPb) 
which represents the envelope functions at A and B sites. 
The eigenstates are labeled by (s,k) with s = +1 and 
— 1 being the conduction and valence bands, respectively, 
and k being the wavevector. The eigenenergy is given 
by £ s k = shvk, and the corresponding wavefunction is 
^sk(r) = e lk ' r F sk /VS with F sk = (e ie ,s)/V2, where k 
and 6 arc defined by (k x ,k y ) = fc(cos#,sin#) and S is 
the system area. 

In a magnetic field B(r) — [V x A(r)] z , the Hamilto- 
nian becomes 7i = Tto + STi. with 8TL = (ev/c)cr ■ A(r). 
The local current density at ro is calculated as the 
expectation value of current-density operator j(i"o) = 
ev cr 5(r — Yq) over the occupied states. In the first order 
perturbation in 8TI, we have 



g v g s e 2 v 2 1 /(e 8 k) - /(e^'k+q) 
x(Fl k <T„F slk+q )(F» a.F.k), (6) 



where g v = g s = 2 are the valley (K,K') and spin degen- 
eracy, respectively, and f(e) = [1 + cxp/3(e — £f)] _1 with 
/3 = l/(fcsT) is the Fermi distribution function. 

At the zero temperature, we can explicitly calculate 
this to obtain 



x(q;£F) 



g v g s e 2 v 1 
16hc 2 g 



0(q - 2k F ) 




2 _■, 2k F 

— sin 

7T q 



(7) 



where kp — \£f\I(Yiv) is the Fermi wave number and 
0(x) is defined by 9(x) — 1 (x > 0) and (x < 0). Signif- 
icantly, x vanishes in range q < 2kp, i.e., no electric cur- 
rent is induced when the external field is smooth enough 
compared to the Fermi wavelength. At ef = 0, particu- 
larly, we have 



x(q;£F = 0) 



g v g s e 2 v 1 
I6hc 2 q 



(8) 



The susceptibility of the carbon nanotube to a uniform 
field perpendicular to the axis has the equivalent expres- 
sion of Eq. ([8]) where o is replaced by 2tt/L with tube cir- 
cumference L. Figure Q] (a) shows a plot of x(q) 
of Eq. ([7]). The susceptibility suddenly starts from zero 
at q = 2k p, and rapidly approaches the universal curve 




FIG. 1: Magnetic susceptibility x(q) in ( a ) graphene and (b) 
conventional 2D system. 



([8]). As a function of ep at fixed q, it is nonzero only in 
a finite region satisfying \ep\ < hvq/2, and its integral 
over ep becomes constant —g v g s e 2 v 2 /(6irc 2 ). Thus, in 
the limit of q — > it goes to 



X(q = 0;e f ) 



g v g s e 2 v 2 
6nc 2 



S(e F ). 



(9) 



This agrees with the susceptibility against uniform mag- 
netic field. 
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Let us consider an undoped graphene (ef — 0) un- 
der a sinusoidal field B(r) — Bocosqx. With the sus- 
ceptibility of Eq. (O, the response current is calculated 
as j(r) = — \g v g s e 2 vB[) / (lQhc)] e y singx. The current in- 
duces a counter magnetic field which reduces the original 
field. The z component of the induced field on graphene 
becomes 



Bind(r) = -a g B(r), a g 



2itg v g s 



16hc 2 



(10) 



Because the ratio is independent of q, Eq. ([TO)) is actually 
valid for any external field B(r), i.e., the magnetic field on 
the graphene is always reduced by the same factor 1 — a g . 
This property holds whenever x(q) is proportional to l/q. 
With the typical value v w 1 x 10 6 m/s, a g is estimated 
as w 4 x 10 -5 , showing that the counter field is much 
smaller than the original. 

The argument of the magnetic field screening can be 
extended in the three dimensional field distribution. Let 
us suppose a situation when a certain magnetic object 
(permanent magnet or electric current) is located above 
the undoped graphene (z > 0), which produces an exter- 
nal magnetic field B(p) in 3D space p = (x, y, z). Then, 
the followings can easily be concluded: (i) On the other 
side of the graphene (z < 0), the induced field becomes 
— a s B(p), i.e., the external field is screened by the factor 
1 — a g . (ii) On the same side (z > 0), the induced field is 
given by a g R z [B(x, y, —z)], where R z is the vector inver- 
sion with respect to z = 0. Namely, this is equivalent to 
a field of the mirror image of the original object reflected 
with respect to z — 0, and reduced by a g . 
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(a) Graphene 



10 r/d 



FIG. 2: Electric current je(r) on graphene induced by a 
magnetic charge q m at z = d. 




0.5 1 
d/a 



For examples, we can calculate the diamagnetic elec- 
tric current and the induced field in several specific ge- 
ometries. We first take a situation where a magnetic 
charge (monopole) q m is located at the point p n = 
(0,0, d), (d > 0) above the graphene plane z = 0. 
The magnetic field perpendicular to the layer on the 
graphene is given by B(r) = q m d/(r 2 + d 2 ) 3 ^ 2 with 
r = (x 2 + y 2 ) 1 / 2 . For Ep = 0, induced m is calculated as 
m(r) = (a g q m /2Tr)(r 2 + ri 2 ) -1 / 2 , and the corresponding- 
current density is given by j(r) = —c{dm/dr) eg with eg 
being an azimuthal unit vector on the xy plane. We note 
that the integral of m(r) over the plane is infinite, show- 
ing that it never looks like a single magnetic dipole even 
if observed from far away. Indeed, the counter field in 
region z > induced by j is the monopole field given by 
a g q m located at — p , as expected from the general ar- 
gument above. The monopole is thus repelled away from 
the graphene with a force of a g q 2 n /(2d) 2 . 

For doped graphenes £f ^ 0, we numerically calculate 
the response current. Figure [2] shows jg{r) for several 
values of kp. We observe the Friedel-type oscillation with 
wavenumber~ 2kp. The overall amplitude exponentially 
decays in region kp > 1/d because the magnetic field 
distribution on graphene has a typical length scale of the 
order of d, while x(<3) vanishes in the long- wavelength 
region such as q < 2k p. 

As another example, we consider line current / par- 
allel to the graphene, which flows along the +y direc- 
tion, passing through the point (0, 0, d) (d > 0) above 
graphene. The z component of the magnetic field on 
graphene is given by B(r) = —2[Iq1x/\c(x 2 + d 2 )]. For 
Ep = 0, the induced current density becomes j(r) = 
— (a g Id/Tr)(x 2 + d 2 ) ~~ 1 e y . The integral of j y in x exactly 
becomes —a g I, i.e., the external electric current induces 
an opposite current on graphene with the amplitude re- 
duced by a g . The magnetic field induced by j in the up- 
per half space [z > 0) becomes equivalent to a field made 
by current —a g I flowing at z = —d, so that the original 
current is repelled by a force a g I 2 /(c 2 d) per unit length. 
When ep 7^ 0, the response current damps for kpd > 1, 
similarly to the case of a magnetic monopole. 

To estimate the diamagnetic force in a possible real- 



FIG. 3: Diamagnetic force per unit area of a semi-infinite 
magnet cylinder with radius a as a function of the distance d 
from the tip to (a) graphene and (b) conventional 2D system. 

istic situation, we consider a case where a semi-infinite 
magnet cylinder with radius a, having a flat end with 
surface magnetic charge density er m , is placed vertically 
above graphene. In real experiments, a may range from 
nanoscale (~ lOnm) up to macroscopic length scale. Fig- 
ure [3] (a) shows the repulsive force per unit area on the 
magnet surface, as a function of the distance d from the 
graphene with several values of kp. The geometry is il- 
lustrated in the inset. For kp — 0, the force is equivalent 
to that made by its mirror charge a g a m under graphene. 
When d ^> a, it approaches the dotted curve given by 
oc 1/d 2 , the force when the surface charge is replaced 
with a point charge q m — 7ra 2 a m . As d goes down to the 
order of a, it deviates from 1/d 2 and reaches 27ra g <7 2 a at 
d = 0, which is exactly the force between a sheet with 
the magnetic charge density a m and another sheet with 
ot g o m , with an infinitesimal gap. For a m which amounts 
to the surface flux of IT (e.g., neodymium magnet), [3l| 
the force is estimated as 0.16 gram force/cm 2 , which is 
surprisingly large as a force generated by a film only one 
atom thick. 

When kp shifts from zero, the force becomes smaller. 
The tail is truncated at d ~ 1/kp, i.e., the repulsive force 
is lost when the distance exceeds the order of the Fermi 
wavelength. This is due to the absence of x(q) f° r l° n g 
wavelength q < 2kp argued above. Similarly, the value 
at d = also decays when kp exceeds the order of 1/a, 
because at d = the spacial distribution of B on the 
graphene has a typical length scale a. 

The graphene diamagnetism is in striking contrast to 
that of the conventional 2D system. If we apply the sim- 
ilar argument to Hamiltonian H = (p + eA/c) 2 /(2m*), 
the nonlocal susceptibility corresponding to Eq. ((TJ yields 

, s 9v9 s e 2 [/ 4k F\ 3/2 a , o, \ / 
* (q;£F)= 2l^4l 1 -^J %-2M-lj, 

(11) 

with Ep = h 2 kp/(2m*). The plot is shown in Fig. 
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[T] (b). When q < 2k p, this is constant at \o = 
—g v g s e 2 / (24irm*c 2 ), which agrees with the usual Lan- 
dau diamagnetism. In the region satisfying q > 2kp, \ 
decays approximately in proportion to 1/q 2 . In highly- 
doped systems such that kp is much larger than the typ- 
ical length scale of external field B(r), the induced mag- 
netization just becomes m(r) xoB(r). For the case of a 
magnetic charge at z = d, we have m(r) cx d/(r 2 + d 2 ) 3 ' 2 , 
and the integral of m(r) over the plane is now finite. The 
induced magnetic field at the distance R 3> d is thus 
dipole-type decaying in proportional to ~ l/R 3 , in con- 
trast to the monopole-type field ~ l/R 2 in graphene. 

Figure [3] (b) shows the force per unit area of the cylin- 
drical magnet with radius a placed above the conven- 
tional 2D system. At large distance d 3> a, it rapidly 
decreases as cx l/d 3 like in the case of a point magnetic 
charge. The peak value at d = has a typical amplitude 
Xo (J m/ a w hh a factor ~ logfcpa. Apart from the factor, 
the ratio of the force at d — of the undoped graphene to 
that of the conventional 2D metal is given by 27ra g a/xoj 
which is roughly the ratio of the values of x(<?) a t 1 ~ 1/ci- 
When the effective mass of GaAs m* ~ 0.067mo is ap- 
plied to xoi the ratio becomes a/(0.01nm). Thus, in a re- 
alistic dimension, the diamagnetic force of the graphene 
is incomparably larger than that of the conventional 2D 
system. Note that, in doped graphene, the diamagnetism 
disappears when a becomes larger than kp . It should 
be noted that x(°) °f graphene does not approach that 
of conventional 2D even in the high kp limit, because of 
difference between linear and quadratic dispersions. 

The singular diamagnetism of graphene is influenced 
by temperature. At zero doping, we expect that x{0) 
deviates from 1/q in Hvq < fc^T, and the divergence at 
q — is rounded to a finite value oc X/QsbT). In other 
words, the temperature affects diamagnetism when the 
typical length scale exceeds 2nHv /k B T ', which is 50 jum 
at T = IK. We expect that the disorder potential gives a 
rnuujilv similar eli.eets to temperature, where 1 I ie_ energy 



scale of level broadening works as finite ksT. [21 



The magnetism is also contributed by electron spins. 
The spin susceptibility x spm (°) is given by the usual 
density-density response function. ^3, HH At Ef = 0, 
this becomes g v fi 2 B q/ {\Ghv) with the Bohr magneton fj,g, 
in contrast to the orbital susceptibility x° rb oc 1/q in Eq. 
ID). The ratio x 8pin /x° rb ~(«x 0.04nm) 2 , is quite small 
in realistic length scales. We also mention that strongly 
disordered graphite and graphene [HI exhibit ferro- 
magnetism due to spin ordering at lattice defects. [HI 

The diamagnetism can be enhanced by stacking 
graphene films. If we have randomly stacked graphenes 



repulsive force given by an external magnet may give rise 
to a mechanical deformation on the graphene sheet. The 
detailed study of this is left for a future work. 
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where the interlayer hopping is neglected, [30j, [3y, [37[ the 
magnetic field would decay exponentially oc (1 — ct g ) 
with the layer number N. We would have an almost per- 
fect magnetic shield when Na g exceeds 1, which amounts 
to the thickness of 10 /xm with the graphite interlayer 
spacing 0.334 nm assumed. We also expect that a strong 
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